
*******************************************************
*Sample characteristics, effectiveness of randomization
*******************************************************

*Table 1 (metric variables)
tabstat age s1_rec s2_rec, s(mean, sd, n) by(group)

*age
robvar age, by(group)
ttest age, by(group)

*smartphone skills
robvar s1_rec, by(group)
ttest s1_rec, by(group)

*internet usage
robvar s2_rec, by(group)
ttest s2_rec, by(group) unequal


*Table 1 (categorical variables)
tabulate gndr_dummy group, column chi2

tabulate edu2 group, column chi2


********************************************************************************************
*RQ 1: Directional Z-tests, Mann-Whitney U-tests, logistic and negative binomial regressions
********************************************************************************************

*H1: Probe nonresponse (Table 2)

tab group OPQ1_nrsp, row
prtest OPQ1_nrsp, by(group)

tab group OPQ2_nrsp, row
prtest OPQ2_nrsp, by(group)


*H1: Logistic regressions (Table E1)

ssc install fitstat, replace
logit OPQ1_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
fitstat

logit OPQ2_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
fitstat


*H1: Logistic regressions (Average marginal effects -- For interested readers; not shown in paper)

logit OPQ1_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)

logit OPQ2_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)


*H2: Not clearly interpretable answers (Table 2)

tab group OPQ1_notint_without_nrsp, row
prtest OPQ1_notint_without_nrsp, by(group)

tab group OPQ2_notint_without_nrsp, row
prtest OPQ2_notint_without_nrsp, by(group)


*H2: Logistic regressions (Table E2)

logit OPQ1_notint_without_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
fitstat

logit OPQ2_notint_without_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
fitstat


*H2: Logistic regressions (Average marginal effects -- For interested readers; not shown in paper)

logit OPQ1_notint_without_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)

logit OPQ2_notint_without_nrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)


*H3: Word count (Table 2)

tabstat OPQ1_wrdcount_withoutsoftnrsp, s(mean, sd, median, n) by(group)
ranksum OPQ1_wrdcount_withoutsoftnrsp, by(group)

tabstat OPQ2_wrdcount_withoutsoftnrsp, s(mean, sd, median, n) by(group)
ranksum OPQ2_wrdcount_withoutsoftnrsp, by(group)


*H3: Negative binomial regressions (Table E3)

nbreg OPQ1_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
nbreg, irr
estat ic

nbreg OPQ2_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
nbreg, irr
estat ic


*H3: Negative binomial regressions (Average marginal effects -- For interested readers; not shown in paper)

nbreg OPQ1_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)

nbreg OPQ2_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)


***************************************************************
*RQ 2: Mann-Whitney U-tests, zero-truncated Poisson regressions
***************************************************************

*H4: Number of themes mentioned 

tabstat OPQ1_No_themes, s(mean, sd, median, n) by(group)
ranksum OPQ1_No_themes, by(group)

tabstat OPQ2_No_themes, s(mean, sd, median, n) by(group)
ranksum OPQ2_No_themes, by(group)


*H4: Zero-truncated Poisson regressions (Table E4)

tpoisson OPQ1_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high, ll (0) vce(robust)
estat ic

tpoisson OPQ2_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high, ll (0) vce(robust)
estat ic


*H4: Zero-truncated Poisson regressions (Average marginal effects -- For interested readers; not shown in paper)

tpoisson OPQ1_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high, ll (0) vce(robust)
margins, dydx(*)

tpoisson OPQ2_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high, ll (0) vce(robust)
margins, dydx(*)


***************************************************************
*RQ 3: Mann-Whitney U-tests, ordinal logistic regressions
***************************************************************

*H6: Survey evaluation interest (Table 5)

tabstat interest_s7rec, s(mean, sd, median, n) by(group)
ranksum interest_s7rec, by(group)


*H6: Survey evaluation difficulty (Table 5)

tabstat difficulty_s8, s(mean, sd, median, n) by(group)
ranksum difficulty_s8, by(group)


*H6: Survey evaluation length (Table 5)

tabstat length_s9rec, s(mean, sd, median, n) by(group)
ranksum length_s9rec, by(group)


*H6. Evaluation length without complete nonresponse (both OPQ1 and OPQ2 not answered)

tabstat length_s9rec if OPQ1_nrsp_hard == 0 & OPQ2_nrsp_hard == 0, s(mean, sd, median, n) by(group)
ranksum length_s9rec if OPQ1_nrsp_hard == 0 & OPQ2_nrsp_hard == 0, by(group)


*H6: Ordinal logistic regressions (Table E5)

ologit interest_s7rec group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high

ologit difficulty_s8 group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high

ologit length_s9rec group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high


*H6: Ordinal logistic regressions (Average marginal effects -- For interested readers; not shown in paper)

ologit interest_s7rec group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)

ologit difficulty_s8 group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)

ologit length_s9rec group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
margins, dydx(*)


*****************************************************************************************
*Alternative regression models for H3 and H4 - For interested readers; not shown in paper
*****************************************************************************************

*H3: Poisson regressions

poisson OPQ1_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
poisson, irr
estat ic

poisson OPQ2_wrdcount_withoutsoftnrsp group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
poisson, irr
estat ic


*H4: Poisson regressions

poisson OPQ1_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
poisson, irr
estat ic

poisson OPQ2_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
poisson, irr
estat ic


*H4: Negative binomial regressions

nbreg OPQ1_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
nbreg, irr
estat ic

nbreg OPQ2_No_themes group_dummy gndr_dummy age edu_dummy_mid edu_dummy_high
nbreg, irr
estat ic


